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ABSTRACT 

We present a generalized version of the Path Interchange Zel'dovich Approximation 
(PIZA), for use with realistic galaxy redshift surveys. PIZA is a linear particle-based 
Lagrangian method which uses the present day positions of galaxies to reconstruct 
both the initial density field and the present day peculiar velocity field. We generalize 
the method by mapping galaxy positions from redshift-space to real-space in the 
Local Group frame, and take the selection function into account by minimizing the 
mass weighted action. We have applied our new method to mock galaxy catalogues, 
and find it offers an improvement in the accuracy of reconstructions over linear theory. 
We have applied our method to the all-sky IRAS Point Source Catalogue Redshift 
Survey. Applying PIZA to PSCz, we are able to obtain the peculiar velocity field, the 
dipole, bulk velocities, and the distortion parameter. From the dipole, we find that 
= 0.51 ± 0.14. We compare the PSCz bulk velocity with that of Mk III and conclude 
that /3 = 0.5 ± 0.15. We compare the PSCz dipole with that of SFI and find that 
= 0.55 ±0.1. 
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1 INTRODUCTION 

Cosmological density and velocity fields may be recon- 
structed using the radial velocity information from a galaxy 
redshift survey. To do this, assumptions must be made about 
the way the peculiar velocity field is produced, the relations 
between the velocity and density fields, and the relation- 
ship between galaxies and the underlying mass distribution. 
The reconstructed fields may then be compared with the di- 
rectly measured fields, e.g. from a peculiar velocity survey, 
and used to place constraints on cosmological parameters. 

The structures seen in the Universe today formed by the 
growth via gravitational instability of small density pertur- 
bations present in the early Universe. The growth and evolu- 
tion of these structures may be studied using methods of two 
types: Eulerian methods that solve the gravitational insta- 
bility equations (see Peebles 1980); and Lagrangian methods 
that follow individual galaxy displacements. These methods 
may be used to deduce the velocity field from the density 
field and vice versa. 

Eulerian methods include those which use linear pertur- 
bation theory applied to the gravitational instability equa- 
tions. The iterative reconstruction methods of Yahil et al. 
(1991), who studied the 1.2Jy redshift survey, and K aiser 
ct al. (1991) and Taylor and Rowan- Robinson (1993) who 



studied the QDOT redshift survey, use this method to re- 
construct the peculiar veloci ty field from galaxy redshift sur- 
veys. Branchini et al. (1999) use a similar method to recon- 
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struct the peculiar velocity of the PSCz. 

Lagrangian methods include the Zel'dovich approxima- 
tion (ZA; Zel'dovich 1970), which was suggested as a way 
to study the evolution of structure in the nonlinear regime. 
Within the ZA the linear displacement field of galaxies is 
extrapolated from initial, Lagrangian coordinates to final, 
Eulerian coordinates. The galaxy peculiar velocity is simply 
the time derivative of the linear displacement. 

The use of Hamilton's principle, that the action is min- 
imized during the evolution of a collection of particles under 
gravity, was introdu ced a s a reconstruction method to cos- 
mology by Peebles ( 1989) ) and developed by Peebles (1994, 
1995) , Shaya, Peebles and Tully ( |l995|), Giavalisco et al. 
( 1993 ) and Susperregi and Binney ( |199^ ) . Sharpe et al. 
( 199S ) recently used Peebles' numerical action method to 
reconstruct the peculiar velocities of PSCz gala xies w ithin 
cz = 2000kms _1 , while Monaco and Efstathiou dl99S|) have 
developed ZTRACE, an iterative reconstruction method 
based on Lagrangian perturbation theory. 

A method combining the ZA and the least action prin- 
ciple is the Path Interchange Z el'doy ich Approximation 
(PIZA) of Croft and Gaztanaga fll997| ) (hereafter CG97). 
They showed that using the ZA the least action solution 
is that which minimizes the total mean square particle dis- 
placement between initial and final galaxy positions. Hence, 
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given a set of final galaxy positions and a random distribu- 
tion of initial coordinates, the least action trajectories are 
those that minimize the total displacement. 

Previously PIZA has been used to reconstruct the ini- 
tial density fields from simulated real space density fields 
(CG97; Narayanan and Croft 1999). Here we wish to apply 
PIZA to a redshift survey in order to reconstruct the pecu- 
liar velocity and real space density fields. We apply our new 
PIZA method to the Point Source Catalogue Redshift sur- 
vey (PSCz; Saunders at al 2000), a survey of 15,000 IRAS 
galaxies to a flux limit of 0.6 Jy at 60pm, covering 84 per 
cent of the sky. 

The outline of the paper is as follows. We review PIZA 
in Section ^. In Section]^ we describe the problems encoun- 
tered when applying PIZA to redshift surveys, and the gen- 
eralizations we have made to the method in order to over- 
come these problems. In Section^ we test PIZA using PSCz- 
like simulations, and compare the predictions for the fields 
with those given by linear theory. In Section 5j we apply our 
new method to the PSCz survey. In Section p we discuss our 
results and present our conclusions. 



2 PATH INTERCHANGE ZEL'DOVICH 
APPROXIMATION 

Here we briefly review the original PIZA method. For m ore 
details, we refer the reader to Croft & Gaztanaga (1997). 



We define a galaxy position in terms of its initial posi- 
tion and displacement: 



r(q,t) =q + £(q,t), 



(1) 



where r is the Eulerian coordinate, q is the initial La- 
grangian coordinate and £ is the galaxy displacement. 

CG97 (§2.1) showed that the action is proportional to 
the mass-weighted squared particle displacement, 



(2) 



where mi is the mass of the i particle. The least action 
principle then requires the minimization of this quantity. 
In CG97 it is assumed that the galaxies are of equal mass, 
simplifying the minimization of equation (|^) . Assuming the 
final, Eulerian galaxy positions are known, the initial posi- 
tions are drawn from a homogeneous distribution and as- 
signed at random to the Eulerian galaxy positions. We shall 
refer to the initial objects as PIZA particles. It is not neces- 
sary that there is only one particle per galaxy; CG97 found 
that increasing the number of particles per galaxy improved 
the accuracy of the reconstruction. Applying PIZA is very 
simple. Two galaxies are then picked at random and their 
PIZA particles interchanged. If this decreases S the inter- 
change is kept, and if not the particles are swapped back. 
Interchanges are attempted until the rate of change of S 
with attempted swaps, the cooling rate, is very low. 



3 PIZA WITH REALISTIC GALAXY 
REDSHIFT SURVEYS 

The PIZA method cannot be immediately applied to a re- 
alistic galaxy redshift survey. In the first instance the radial 



positions of galaxies are estimated from their redshift, which 
introduces distortions in galaxy positions due to the pecu- 
liar velocities of galaxies. This distortion will depend on the 
rest frame of the observer, which we shall assume is the Lo- 
cal Group frame. To apply PIZA we must first deal with 
this distortion. Secondly most galaxy redshift surveys are 
flux limited, so that the observed number of galaxies drops 
off with distance as only the brightest galaxies can be seen. 
This affects the number density of galaxies, which must now 
be regarded as sampling the total galaxy population. This 
means that the galaxies must be assigned effective masses 
which must be included in the PIZA minimization. Finally, 
the sky coverage of galaxies may not be complete and some 
way must be included in the reconstruction to account for 
regions of incompleteness. In this Section we shall discuss 
each of these effects in turn. 



3.1 PIZA in redshift space 

The peculiar motions of galaxies, v, mean that in the CMB 
frame the radial coordinate in redshift space is a mixture of 
spatial position and projected velocity, 



s(r,t) = r + (f.v(r,t))f, 



(3) 



where s is the redshift position. In the ZA the velocity field 
of galaxies is a linear extrapolation, and related to the dis- 
placement field by 

v=t = m)t (4) 

where / = rfl n<5/rflna ~ Q , os is the growth rate of linear 
perturbations (Peebles 198C). If galaxies are a locally biased 
representation of the density field then 5 ga i = f[5 maS s] ~ 
b5 mass to first order. We can then define a galaxy displace- 
ment field which is related to the velocity field by 



Pi 



where 



P 



m) 



(•») 



(6) 



is the redshift distortion parameter. The redshift space dis- 
pla cement field is related to t he real space displacement field 
by (|Taylor fc Hamilton 199(1 



G = "P«ti, (7) 

where V%j = 8^ + (3fifj is the redshift space projection 
tensor, and SE is the Kronecker tensor. Inverting equation 



(@)> S ives ( [Taylor 



Valentine 1999| ): 



where 



P 



(8) 



(9) 



is the inverse redshift space projection tensor. Squaring this 
gives us the expression for the linear displacement field in 
terms of linear redshift displacements and (3; 

1 



£ 2 = £ 



(i 



(10) 



(1 + /3) 2 ' 

where = £.r is the displacement vector in redshift space 
projected along the line of sight. Thus we must assume an 
a priori value for f3. 
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3.2 PIZA with Local Group frame redshifts 



Redshift space distortions should also take into account the 
motion of the observer, in this case the Local Group motion. 
The transformation f rom real to redshift space in the Local 
Group rest frame is (Taylor & Valentine 1999) 



Radial Velocity (kms ) 



Radial Velocity (kms ') 



(ii) 



where the super- or sub-script LG denotes a variable mea- 
sured in the Local Group frame, and £ (0) is the displacement 
of the Local Group. The real-space displacement field can 
now be expressed in terms of redshift displacements in the 
Local Group frame as 



G96-(o) + E 



(12) 



(1+/3) 2 

To apply equation (^) we require the displacement of the 
Local Group. We find this by placing a galaxy representing 
the Local Group at the origin and using its displacement as 
the Local Group displacement. 



3.3 Incomplete Sky Coverage 

The incomplete sky coverage may cause problems with 
PIZA, as with other reconstruction methods, since in these 
regions we do not know the true galaxy distribution. A num- 
ber of choices for filling in these regions present themselves. 
One could populate the incomplete regions randomly on the 
sky, assuming no information or leave the regions empty of 
both particles and galaxies. These are equivalent except for 
particles adjacent to the mask. Alternatively, one could in- 
terpolate across using correlations in the density field to fill 
in the gaps. In reconstructing PSCz we use t he Fourier In- 
terpolation Scheme of Saunders et al. (1999) which allows 



optimal, nonlinear interpolation across these regions. 



3.4 PIZA with a selection function 

In a flux limited survey, the number density of survey ob- 
jects decreases with increasing redshift, as galaxies at greater 
distances must be brighter to be seen above the flux limit. 
This is quantified by the survey selection function, <j>(r), the 
expected number of galaxies seen at distance r above the 
flux limit in the absence of clustering. 

In order to account for this change in the number den- 
sity with distance, we pick initial PIZA particles with the 
same selection function as the galaxies. To take account of 
those galaxies that fall below the flux limit we assign each 
galaxy a mass of l/<f>(r) and each PIZA particle a mass of 



m, = l/(u(j>(ri)), 



(13) 



where v is the mean number of particles assigned per galaxy. 
In this work we use 10 particles per galaxy. We use u = 
10 in order to reduce the shot noise in the reconstruction. 
Equations [l^ and |l3| are used in equation ^ to minimize the 
displacement squared. 

With the inclusion of particle masses it is possible for 
S to decrease but light particles to acquire unrealistically 
large trajectories. We have therefore added an additional 
constraint that the rms displacement squared in redshift 
space should also decrease. 
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Figure 1. Typical simulation results: comparison of PIZA and 
linear theory reconstructions. The left hand panel shows PIZA 
radial velocities compared with simulation, the right hand shows 
linear theory compared with simulation. Velocities were smoothed 
with a Gaussian of radius 20/i~ 1 Mpc. The dotted line indicates 
the least-squares fit. 



3.5 Second Order Effects 

The masses we have assigned to the galaxies are based on 
the selection function at their redshift distance. Because of 
peculiar velocities, these can be in error. We therefore de- 
veloped a modified version of PIZA, MASSPIZA, whereby 
after PIZA is run, particles are randomly reassigned from 
galaxies whose mass has been overestimated to those where 
it has been underestimated. PIZA is then rerun and the 
whole process iterated until the masses are acceptable and 
the trajectories minimised. Unfortunately, this procedure is 
in general unstable to the rocket effects described by Kaiser 
(1987). In the section below, we show that PIZA without 
mass reassignment gives acceptable results in comparison 
with simulations; the agreement is in general not improved 
by the second, MASSPIZA stage. Therefore, we have not 
incorporated it in any of the results quoted in this paper. 



4 TESTS ON SIMULATIONS 

In order to test the new PIZA method, we have applied it 
to a set of PSCz-like mock catalogues. For more details of 
thes e mo ck catalogues, we refer the reader to Branchini et 
al. (200C). The catalogues are of two CDM cosmologies, a 
Sl m = 1, P = 1, with T = 0.25 and a flat 
= 0.3 and SIa ~ 0.7, = 0.50 and are lim- 
We test the reconstruction of the 



critical model 
model with Q„ 
ited to v < 15000fcms" 
radial peculiar velocity field, comparing it with that of the 
simulation. To compare the fields the peculiar velocities were 
binned in real space and smoothed with a Gaussian smooth- 
ing kernel. We have also compared the PIZA reconstructions 
with the linear theory reconstructions for the mock cata- 
logues from Branchini et al. Figure [j] shows typical results 
comparing PIZA and linear theory reconstructions with the 
simulations. The dotted line indicates the the least squares 
fit to the data. As stated above, we must assume an a priori 
value for /3. We find that the reconstruction works well for 
any sensible value of (3. We find that generally the gradient 
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Scatter 
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0.76 


0.10 
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0.95 


0.27 
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0.61 


0.14 


A0.7 4tt 


0.43 


0.05 


14 


0.74 


0.04 


T0.25 4?r 


0.84 


0.18 


13 


0.65 


0.13 



Real Space 



Table 1. Results from tests on simulations. AO. 7 denotes the Qa 
model, and TO. 25 the critical model. P denotes simulations with 
PSCz sky coverage and 47r simulations with full sky coverage 



reconstructed by PIZA is less biased than that reconstructed 
by linear theory but that the scatter is somewhat worse. This 
is probably one of the limitations of PIZA as the scatter in 
redshift space reconstructions produced by PIZA is quite 
large (see also CG97). 

We find that PIZA works better for low density models 
as they are less dynamically evolved. 

The accuracy of the reconstruction varies between dif- 
ferent sets of initial conditions (i.e. PIZA particles) due to 
the random positioning of particles. There is also a scatter 
between reconstructions using the same initial conditions 
due to the random swaps made by PIZA. This probably re- 
flects that there is not a simple minimum for the action, 
S, but that the configuration space is complicated by the 
particle discreteness. 

We estimate f3 from the trajectory of the Local Group 
galaxy using equation |B| and assuming the magnitude of the 
CMB dipole. Table |l| shows the reconstructed /3, the dipole 
misalignment angle, and the slope of the peculiar velocity 
comparison, with the scatter due to different initial con- 
ditions. We have tested the effects of having the Galactic 
Plane filled or empty. This is important for PSCz as we know 
there are large structures, such as the Great Attractor, in 
the Plane. With the simulations we have tested PIZA using 
An sky coverage, and with the PSCz mask applied. With the 
mask empty we find that the average misalignment angle is 
24° , and with full-sky coverage it is 14° . 



5 RECONSTRUCTING PSCZ 

We have applied our new PIZA method to the PSCz survey 
(Saunders et al. 2000) out to R = 300/i _1 Mpc, using the 
parametrized IRA S selection function of Mann, Saunders 
and Taylor (1996). For the reconstructions we have filled 



the masked regions with interpolated galaxies. 

Figure (3 shows the reconstructed PSCz velocity field in 
a slice 20 h Mpc thick centred on the Supergalactic Plane. 
Figure ^ shows the same slice out to 50 /i _1 Mpc. The arrows 
show the galaxy trajectories, with the arrows starting at 
the initial positions and the arrow heads at the galaxy real 
space positions. Generally the velocity field at the edge of 
the volume is quiet, apart from flow towards Shapley at 
(120,90) and the void at (0,-120). The Coma cluster is seen 
at (0,70). The Cetus wall extends northwards from (0,-100) 
towards Perseus-Pisces at (50,-20). Virgo is seen at (0,10). 
The Great Attractor is at (-40,20), and we see little evidence 
of backfall on to it on the opposite side to the G.A. from the 
Local Group. 
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Figure 2. PSCz velocity field within 150/i _1 Mpc 
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Figure 3. PSCz velocity field within 50ft _1 Mpc 
5.1 Cosmological dipole 

We find an average dipole direction of (l,b)=(264.4°,41.7°). 
We find that j3 = 0.51 ± 0.14. This error has 3 sources: ran- 
dom error due to different initial conditions; misalignment 
between the reconstructed and CMB dipoles; cos mic vari- 
ance and shot noise on the reconstructed dipole (Taylor &i 



Valentine 1999). A similar result is found by only using the 



z component of the dipole. 



5.2 Radial velocity field 

The reconstructed dipole and bulk flows are shown in Fig- 
ure ^ We have compared our bulk flows with those from 
Mark III and the linear theory PSCz reconstruction (Bran- 
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Figure 4. PSCz bulk velocities and dipole within 150/i _1 Mpc 
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Table 2. SFI vs PSCz peculiar velocities with respect to shells 
of galaxies. PSCz reconstruction assumes /3 = 0.5. 



chini et al. 1999), and conclude that (5 = 0.5 ± 0.15. This is 
consistent with the results of Tadros et al (1999) based on 
the measurement of the degree of redshift space distortion 
in the PSCz. 



5.3 Comparison with SFI 

We have undertaken a comparison of the reconstructed 
PSCz velocity field with the SFI data set of peculiar ve- 
locities of spiral galaxies (Giovanelli et al. 1998). We have 
calculated the average peculiar velocity, with respect to the 
Local Group, of galaxies in redshift slices, to be as closely 
as possible comparable to Giovanelli et al.. We find reason- 
able agreement in amplitude, though there are significant 
differences in direction, for (3 = 0.55 ± 0.1. Table H shows 
the average peculiar velocity in shells for PIZA and SFI. 



6 DISCUSSION 

We have presented a generalized version of the PIZA method 
for reconstructing velocity and density fields from galaxy 
surveys. These modifications to the original method are nec- 
essary if it is to be applied to realistic galaxy redshift surveys 
with a selection function. We have used the inverse redshift 



space operator of Taylor and Valentine (1999) to write the 
mean square particle displacement in terms of the redshift 
space positions of galaxies, and to map redshifts from the 
Local Group rest frame to real space. The PIZA particles 
obey the same selection function as the galaxies. Doing this 
we can reconstruct the real space density field, the galaxy 
initial positions and the peculiar velocity field. 



We have tested the method using PSCz-like CDM simu- 
lations, testing the reconstruction of the real space positions 
and the peculiar velocities. We have also compared our re- 
constructions with linear theory reconstructions of the sim- 
ulations. We have applied our method to the PSCz survey, 
and reconstructed the local velocity field out to 300 /i _1 Mpc 
from the Local Group, and have calculated the bulk flow and 
the dipole. We find the reconstructed dipole to be 16° ± 4 
away from the CMB dipole. From comparison between the 
PIZA reconstruction and the Mark III bulk flow, we find 
that (3 — 0.5 ± 0.15, while from a dipole analysis of the 
PIZA reconstruction, we find that (3 = 0.51 ± 0.14. 
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